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Abstract 

In this paper we propose a heuristic technique for distributing points 
on the surface of a unit n-dimensional Euclidean sphere, generated as 
the orbit of a finite cyclic subgroup of orthogonal matrices, the so called 
cyclic group codes. Massive numerical experiments were done and many 
new cyclic group codes have been obtained in several dimensions at var¬ 
ious rate. The obtained results assure that the heuristic approach have 
performance comparable to a brute-force search technique with the advan¬ 
tage of having low complexity, allowing for designing codes with a large 
number of points in higher dimensions. 
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1 Introduction 

Consider a unit vector x in the n-dimensional Euclidean space K” and a set of M 
distinct orthogonal n x n matrices G = {51,52, ’ , 5 m }- The orbit of x under 

the action of G is the set of points C = {giX,y gi G G}. Due to orthogonality, 
the points in C belongs to the surface of the unit sphere 5”“^ C M”. In the 
case of the set G is a finite cyclic multiplicative subgroup, the orbit C is called 
a cyclic group code. 

Cyclic group codes belongs to the family of Slepian Group Codes, introduced 
in [1] or, in a more general sense, geometrically uniform codes [2]. Such codes 
have been widely applied in communication theory, for instance to match signal 
sets to groups [3]. In general, point sets on the unit sphere are useful for 
communicating over a Gaussian channel and are a natural generalization of 
phase shift keyed signal sets (PSK) to dimensions greater than two |4l[5]. 
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Improve the design of group codes in higher dimensions may have a consider¬ 
able impact in error correcting codes and digital communications. In particular, 
it was proved in |6] that the Shannon capacity of certain important channels, 
as the AWGN channel with m-PSK modulation, can be achieved using commu¬ 
tative group codes. As we show in Section]^ cyclic group codes may asymptot¬ 
ically approach the sphere packing bound for commutative group codes, while 
have a simpler structure. For instance, all the M matrices in a cyclic subgroup 
of n X n orthogonal matrices can be generated from only one [n/2J-dimensional 
integer vector (see Q in Section]^. 

For any fixed subgroup G, the initial vector problem (IVP) is to choose 
X which maximizes the Euclidean distance between any two points on the orbit 
C. This problem has not been solved for general groups, but rather it is known 
only for special cases such as cyclic groups [7], commutative groups [S], finite 
reflection groupsj^ and permutation codes [TU]. We remark that, for the cyclic 
and commutative cases, the IVP can be stated as a linear programming problem 
as reviewed in Section I01 

An important quota in the theory of spherical codes is the maximum number 
of points displayed on the surface of a unit sphere with distance not smaller 
than a certain d. In El it is derived the bound Q on the maximal number 
M of points in a (2fc)-dimensional commutative group code as a function of 
the minimal distance d and the maximal center density of a lattice packing in 
fc-dimension. In this case, 

M <-— * < -^— fl) 

(arcsin^)'' (arcsin ’ 

where Ak is the maximal center density of a lattice packing in and x = 
{xi,X 2 , • ■ ■, X 2 n) is the initial vector. 

Despite the relevant literature on group codes (indeed spherical codes in 
generaQ, explicit construction of cyclic group codes for a large number of points 
is still a challenge problem, particularly in high dimensions. It is worth to remark 
that the classical Simplex code is a cyclic group code with n + 1 points in n- 
dimension and the Biorthogonal code is also a cyclic group code with M = 2n 
points in odd dimension {Biorthogonal is not cyclic in even dimension). A more 

thorough survey of spherical codes can be found in m and explicit constructions of 
spherical codes were proposed in [121ll3lfT4) . Lists of good spherical codes and also packings, 
coverings and designs can be found online at m- 
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detailed discussion can be found in Chapter 8 of m- 

In |S] it is presented a method for finding an optimum (2fc)-dimensional 
commutative group code of a given order M, which is based on the structure 
of lattices related to these codes. Although that method provides a significant 
reduction in the number of non-isometric cases to be analyzed, it still demands 
a brute-force search on about 

M" 

THp{M) 

cases, where stands for the Euler phi function of M. 

It is easy to see that, when M and n increase, the number ([^ become 
computationally prohibitive. For instance, for {2k, M) = (16,1024), which cor¬ 
responds to a commutative group code in C with M = 1024 points, 
the approach proposed in |S] to find the optimal code demands a full search 
in a set with about 2®® elements. In addition, it also means solving 2®® linear 
programming problems to find the best initial vector for each case. 

In this paper, we present a heuristic technique inspired on [S] , for designing 
(2A:)-dimensional cyclic group codes for a given number of points M, avoiding 
the brute-force search and hence allowing for find codes with a large number 
of points in high dimensions. Using this technique we were able to design new 
cyclic group codes for M up to 2^® = 524288 points in dimensions up to n = 48, 
with good performance (Figure]^. 

Although the heuristic can be adjusted for design commutative group codes, 
we restrict our attention in the cyclic case because these codes have a simpler 
and stronger structure and, moreover, the codes founded with our heuristic 
suggests a conjecture that the cyclic family can asymptotically approach the 
sphere packing bound for commutative group codes (Figure]^. 

The rest of the paper is organized as follows: In Section we review fun¬ 
damental concepts on cyclic group codes and set our notation. In Section we 
describe our heuristic method and present a pseudo-code algorithm for that. In 
Section we show various cyclic group codes designed using the heuristic ap¬ 
proach and some performance comparisons. Finally, in Sectionj^we summarize 
and conclude the paper. 

2 Cyclic group codes 

Let On be the multiplicative group of orthogonal matrices nx n and Gn{M) be 
the set of all order M cyclic subgroups in On- 
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A cyclic group code C on the surface of a n- dimensional euclidean sphere is 
a set of M vectors generated as the orbit of an initial vector x G C K" by 

a given G G i.e. 


C := Gx = {gx, g G G} . 

The minimum distance in C is defined as: 

d := min ||x —j/||= min —x||, 

x.y e c Qi e G 

gi=A In 

where 11.11 denote the standard Euclidean norm and In is the identity matrix of 
order n. 

We use C{M,n,d) to denote a code C in K" with M points and minimum 
distance equal to d. 

The minimum distance of a cyclic group code C, depends on the initial vector 
X and also on the orthogonal matrices on G (in other words, the distance also 
depends on the real representation of the subgroup G). For every finite cyclic 
subgroup G it is possible to find a generator matrix g, such that G = {(5)'’}^o^, 
where = /„ is the identity matrix in M". Thus, C = {( 5 )-'a;}^Q^ and 
therefore 

d= min ||(g)'’x — x||. 

3 7^ 0 

For every commutative group code of orthogonal matrices there is a well 
known real-irreducible representation, which can be stated as follows: 

Theorem 2.1 (JT^ Theorem 12.1) Every commutative group G G can 

be carried by the same real orthogonal transformation q into a pseudo-diagonal 
form: 


<19,(f = [-^*( 1 ), • ■ ■, Ri{k),p.{i)2k+i, • ■ ■, Ti.i)n]r 


where Ri(bj) = 


/ 2'Kbin \ . / 2'Kbi-i \ 

cos(Tr^) -sm(^) 

• / 2'Kbin \ / 27zbi-i \ 

cos(^) 


(3) 


bij G Z, 0 ^ bij ^ M and p.{i)i = ±1, Z = 2/c -|- 1,..., n, j = 1,..., fc, G G. 
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This result holds, in particular, for the generator g, which means that every 
cyclic group code in even dimensions is equivalent to one in which the generator 
matrix g has a quasidiagonal form 


g = [R{bi),Ribn/2)]nxn, where R{bj) 




(4) 


Note that each 2x2 block R{bj) in Q is a rotation matrix by an angle of 
{2TTbj /M). Thus the generator matrix g of any cyclic group G G On= 2 k is defined 
by a vector 

b = ( 6 i, • • • ,bk) with 0 < bi < M and gcd(&i, 62 , • • • , bk,M) = 1 that repre¬ 

sent the rotation blocks. In addition, as pointed out in [S] up to conjugacy we 
can consider 0 < bj < M/2. 

If n is odd, then 

9 ~ • ■ • , )? 9]nxnj 

where /r = I if M is even and ^ = — I if M is odd [ 3 . Thus, the generator g 
may be always defined by the vector b. 


2.1 The Initial Vector Problem for cyclic group codes 

Using the decomposition Q it is possible to state the initial vector problem 
(IVP) for commutative group codes (hence for cyclic) as a linear programming 
problem as follows |318]: 

For a given G G On the initial vector problem is to find x G 5'"“^ such that: 


max 


( minll^iX — X 

V *5^0 


From Theorem 2.1 and Equation (|^ we can define 
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Trh 


fi{x) := \\g^x-x\\ =4> (a: 2 ,_i-kX 2 ,-)sin ^+4 > x., i = ^ 


9=1 


M 


E 

1 = 9+1 


,M. 


Let 


Vj = 




X 


2 

j + fc 


, if j = l,...,fc 
, if j = fc -I- 1,..., n — A: 
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we have that 


My) = 


i=i 


2 

~w 


n—k 

j=/c+l 


and, therefore 


max (min\\g^x - x\M = max (min/j( 2 /)), where ,j/= (yi, •••, j/fc). 
xeS"--! y6S<"-'=-'-) 


This max-min problem can be stated as a linear programming problem as 
in 0: 


min c*w 


where: 


s.t. 


Aw < b 
w >0 



A = 


(-1,0,...,o) e 
(i,-i,o,...,o) G kL“J+2 
/ 0 1 ... 1 \ 

0 -1 ... -1 

1 mi,i ... mn-q,l 


rriij = is the {j — th) coeficient of/^. 


Remark: since this LP-problem has {n — k + 1) variables and [M/2J con¬ 
straints, when M » n it is worth to solve this by using the dual-simplex 
method. This is the method we have used for running our numerical experi¬ 
ments presented in Section]^ 


2.2 On optimal cyclic group codes 

From the previous discussion we can conclude that for a given number of points 
M it is possible to run a brute-force search method for finding an n-dimensional 
cyclic code by analyzing all possible vectors b and solve the correspondent IVP 
for each one of them. In fact this approach was used in [7] for find some optimal 
cyclic group codes for small number of points M and improved for commutative 
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group codes in [S]. 

The main weakness of the brute-force approach is that it leads to check 

/M/2\ 

about cases and, for each one of them, solving the IVP to maximize 

the minimum distance. Indeed, it has been shown in that it is possible to 
reduce the number of cases to less than ([^, by considering isometries arising 
from lattices. Using this reduction, it was presented in [S] several optimum 
commutative group codes in dimensions 4 and 6 for M < 1000. Neverthless, as 
pointed out in the introduction, the number (H may become computationally 
prohibitive when M increases. 

In the next section we describe a heuristic approach that allowed us to design 
codes for n up to 48 with M < 2^® points. As we can see in Table the 
codes found by the proposed heuristic have comparable performance with some 
optimal codes known and are possibly the best published lower bounds for many 
values of (M, n). 


3 A heuristic approach for designing cyclic group 
codes 

We start with the same approach as in |S], i.e. by exploring the connection 
between cyclic group codes and lattices. Formally, let C be a cyclic group code 
in generated by a group G S 0„, G = {(3)^ }^o^i where g is defined as 

in Q by an integer vector 6 = (6i, • • • ,bk)- 
The associated lattice Aq is defined by 

Ag := {(6 i,..., bk) G : [R{bi), ..., R(bk)] G G} , 

where R{bj) denotes the rotation in by an angle of 27Tbj/M and [ii(6 i),..., R{bk)] 
denotes a pseudo-diagonal matrix, according to (@. 

As pointed out in j8 ] Ag contains MZ^ := {M{zi, Z 2 , ■ ■ ■, Zk), Zi G Z} as a 
sub-lattice. Thus, inside the hyperbox [0, M)^ there are exactly M points of 
Ag, which correspond to representatives of the elements of G, i.e., 

[0, M)” D mod M, i = 1, 2,..., A/} . 

If a; = (xi, 0, a; 2 , 0,..., Xk-, 0) is an initial vector for the code C, we can also 
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define a lattice Aa{x) by 


r 1 1 

I 2TrX2 I 

^ 0 ( 2 ;) := < , b:bGAc\. 

il J 

Under this assumption the code C is the image tfjx{Ac{x)) C where 


ipx{y) = ( Xi cos ( ^ j , Xi sin f ^ j ,..., a;fc cos f ^ j , Xk sin f ^ ) ) (5) 


is the standard parametrization of the torus with radii Xi, i = 1,..., k mi. 

As an example, let us consider two 4-dimensional cyclic group codes Ci and 
C 2 , both with M = 25 points, defined by vectors 61 = (1,2) and b = (3,11), 
respectively. These two particular codes are isometrics and so their associated 
lattices Ai and A 2 are isomorphic, as illustrated in 




10 15 20 



Figure 1: The associated 2-dimensional lattices Ai and A 2 of two isometric 4- 
dimensional cyclic group codes C\ and C 2 defined by the vectors bi = (1, 2) and 
b = (3,11), respectively. 


This association between a 2A:-dimensional cyclic group codes and a k- 
dimensional lattices is behind our heuristic method, as we describe in the next 
Section. 






3.1 The heuristic 


The key idea: 

When we analyze the lattice Ai, plotted on the left side of Figure which 
correspond to the group defined by the vector b = (1, 2), we may conjecture that 
the minimal distance in the 4-dimensional group code Ci, is given by ||'0x(0,0) — 
■i/jx (1, 2) 11. Thus, if we want to find codes with large minimal distances, we should 
discard candidates such that || 6 || is to small. On the other hand, we must hold 
our greed, because if || 6 || is bigger than a certain threshold, then the minimal 
distance of the group code will be defined by other pair of points, rather then 
||V^a:(0,0) — ' 02 :( 1 , 2 )||, as occurs, for instance, in the lattice A 2 (right hand side 
of Figure . 

In what follow we derive a fair threshold for || 6 || in order to get good cyclic 
group codes, which is the basis for our heuristic approach. 

Since vector b = {bi,... ,bk) defines a generator g of a cyclic subgroup G, 
if and only if, gcd( 6 i,..., &fc, M) = 1, we can trivially guarantee this condition 
by choosing 61 = 1 and, of course, reducing the set of candidates (eventually 
loosing optimality). So in the heuristic we assume that 61 = 1 and it is the first 
reduction comparing to a force-brute approach. 

The main reduction in the set of candidates comes from writing the inequality 
Q in terms of d, 

( 6 ) 




and get a target (upper bounded) distance d := 4 sin ( ) in the cyclic 

group code C{M,n). 

We know that A^ = ^ ^ where llwminll and vol(A) are the minimal 

2''vol(A) 

norm and the volume of the densest fc—dimensional lattice, respectively. Thus, 
for lattices associated to cyclic groups we can approach ||wmin|| = ||^>|| and 
consider vol(A) = . 

Thus, we may restrict our search on candidates b such that bi = 1 and 


2My/k 


arcsin(d/4) 


(7) 


In Figure we show an example of the choice of candidates according this 
heuristic for finding a 4-dimensional cyclic group code with M = 20 points. The 
circle (arc) has radius d =~ 4.805. In this case, the only vector to be tested is 
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Figure 2: The associated lattices in of a 4-dimensional cyclic group codes 
with M = 20 points defined by the vector b = (1,5) and an arc-circle of radius 
d. 

b = (1,5) which has || 6 || ~ 5.1 ~ d. This vector corresponds to a cyclic group 
code with minimum distance equal to 1.19 while the best commutative group 
code for {M,n) = (20,4) found in [S] has distance equal to 1.24 and the bound 
©is equal to 1.465. 

We remark that, in the 2-dimension, the only choice we have is on coordinate 
62 , but when the dimension increases, it become more flexible, allowing testing 
more candidates b and hence increasing the chance of obtaining good codes. 

Algorithm 1 summarize a pseudo-code for the heuristic and some results are 
described in the next Section. 


4 Results 

In this section we present numerical results obtained with the heuristic described 
in the previous section. Various codes were found in dimensions up to 48 for M 
up to 2 ^®. 

In Figure we present a compilation of several cyclic group codes found 
in this paper. The results have been plotted, as in |7], in a Jabobs’ efficiency 
chart for comparison. In this chart, every code C(M, n) with minimal distance d 
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Algorithm 1: Heuristic method for design cyclic group codes 


inputs : The number of points M; The dimension n = 2k; The maximal 
center density of a lattice packing in A^; A parameter Q to 
define the number of cases to be tested 
outputs: A vector b, corresponding to the generator of a cyclic group 
code C{M,n); The optimum initial vector x and minimum 
distance d. 

begin 

initialization; 

dbest Oj 
^best ^ 0 ; 


d ^ 4 sin 


. / Tr'^Ak 

2M^ 


W < -arcsin(c?/4); 

TT 

while C > 0 do 

1) Chose a random vector v = (vi,..., Vk), such that: 

/ — vf — ■ ■ ■ — vf 

vi = i,Vi ^ Uj+i < y- 

Vk = ^W^-vl - ; 

2) Solve the LP corresponding to the IVP to get x and the 
minimum distance dmin] 

if d^i^fi > dff^gi th-Gii 
dbest d-iYiin ; 

^best 

^best 

_ C^C-1; 

d dbesti 
^ '^best ) 





and minimal angle separation p = 2 arcsin(d/2) is represented by a point (i?, K), 
where 

K = {1 — p) log 2 M and R = ^ log 2 M. This charts exhibit codes in dimen¬ 
sions 8,10, • • • , 24 for M = 2®, 2^, • • • , 2^® points. The dotted and dashed lines 
correspond to the performance of biorthogonal and simplex codes, respectively. 

K 


48-codes 

t 

/ 44-codes 



0.5 1.0 1.5 2.0 2.5 W 


Figure 3: Efficiency comparison for commutative group codes, for 1000 cases 
analysing. 

As we can see, most of the codes presented in outperform the classical 
Simplex and Biorthogonal codes, and also the codes showed in [7]. Therefore, 
using the heuristic we were able to design cyclic group codes in high dimension 
with a large number of points that really have good performance. We believe 
the codes presented in this paper might be the best known lower bounds for 
minimal distances of cyclic group codes in several dimensions published up to 
now. 

In Figurej^we compare the ration between the minimal distance d of various 
cyclic group codes found by using the heuristic and the target distance d = 
dbound obtained from the sphere packing bound In the x-axis we plotted 
log 2 (M) and the y-axis show i- As we can see, when M increases the ratio 
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M 

Q = 512 

Q = 1024 

Q = 1536 

Q = 2048 

64 

1.27269 

1.28386 

1.2859 

1.29514 

128 

1.21733 

1.21962 

1.23481 

1.23481 

256 

1.15941 

1.16528 

1.16905 

1.23721 

512 

1.13819 

1.19334 

1.19334 

1.19334 

1024 

1.11273 

1.12457 

1.12457 

1.12457 

2048 

1.09191 

1.09191 

1.09191 

1.09191 

4096 

1.06506 

1.06506 

1.07339 

1.07339 

8192 

0.99974 

1.02577 

1.02577 

1.03603 

16384 

1.02518 

1.03103 

1.03103 

1.03103 

32768 

0.98699 

1.01643 

1.01643 

1.01643 

65536 

0.99445 

0.99445 

0.99445 

0.99445 

131072 

0.94748 

0.94748 

0.96942 

0.96942 

262144 

0.95448 

0.95448 

0.95565 

0.95966 

524288 

0.91805 

0.95606 

0.95606 

0.95606 


Table 1: Heuristic performance depending on the number of candidates Q 
tested for each case. Examples for designing 48-dimensional codes whit M = 
26,2^,-- - ,219. 


M 

Bound 

Optimum 

Heuristc 

10 

1.474 

1.224 

1.224 

20 

1.054 

0.959 

0.917 

30 

0.864 

0.831 

0.769 

40 

0.750 

0.714 

0.707 

50 

0.672 

0.628 

0.609 

100 

0.476 

0.468 

0.433 

200 

0.337 

0.330 

0.317 

300 

0.275 

0.273 

0.259 

400 

0.238 

0.237 

0.221 

500 

0.213 

0.211 

0.200 

600 

0.194 

0.193 

0.180 

700 

0.180 

0.180 

0.167 

800 

0.168 

0.168 

0.162 

900 

0.159 

0.158 

0.148 

1000 

0.150 

0.149 

0.146 


M 

Bound 

Optimum 

Heuristc 

10 

1.820 

1.414 

1.345 

20 

1.465 

1.240 

1.190 

30 

1.287 

1.133 

1.056 

40 

1.173 

1.044 

1.007 

50 

1.091 

0.976 

0.946 

100 

0.870 

0.804 

0.786 

200 

0.692 

0.673 

0.633 

300 

0.605 

0.585 

0.568 

400 

0.550 

0.540 

0.525 

500 

0.511 

0.504 

0.479 

600 

0.481 

0.472 

0.458 

700 

0.457 

0.445 

0.439 

800 

0.437 

0.427 

0.415 

900 

0.420 

0.413 

0.403 

1000 

0.406 

0.397 

0.394 


Table 2: Comparison between codes found by the heuristic for several values of 
{M, n) with respective optimum cases and bound. 


approaches to 1, but the convergence is slower as the higher is the dimension. 
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logaM 


Figure 4: Comparison between the minimal distance of cyclic group codes found 
by using the heuristic and the sphere packing bound for commutative group 
codes. 

5 Conclusions 

We presented a heuristic technique for designing n = 2fc-dimensional cyclic 
group codes for given number of points M. The heuristic explores the relation 
between such codes and fc-dimensional lattices. Numerical experiments were 
done and many new cyclic group codes have been obtained in several dimensions 
at various rate. The obtained results assure that the heuristic approach have 
performance comparable to a brute-force technique with the advantage of having 
low complexity, which allows for designing codes with large number of points 
in higher dimensions. The results also suggest that the sphere packing bound 
for commutative group codes can be asymptotically approached by cyclic group 
codes. 
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